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Doppler Effect of Nonlinear Waves and Superspirals in Oscillatory Media 
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Nonlinear waves emitted from a moving source are studied. A meandering spiral in a reaction- 
difltusion medium provides an example, where waves originate from a source exhibiting a back-and- 
forth movement in radial direction. The periodic motion of the source induces a Doppler effect that 
causes a modulation in wavelength and amplitude of the waves ("superspiral"). Using the complex 
Ginzburg-Landau equation, we show that waves subject to a convective Eckhaus instability can 
exhibit monotonous growth or decay as well as saturation of these modulations away from the 
source depending on the perturbation frequency. Our findings allow a consistent interpretation of 
recent experimental observations concerning superspirals and their decay to spatio-temporal chaos. 

PACS numbers: 82.40.Ck, 05.45.-a, 47.54.+r 



Introduction. - Periodic nonlinear waves are a trade- 
mark of nonequilibrium systems 0. In one dimension 
(ID), they can appear in systems with periodic bound- 
ary conditions (BCs) as well as in open geometries, where 
BCs select a unique pattern In two dimensions (2D), 
rotating spiral waves are frequently observed. Therein, 
periodic waves emerge from the region of the spiral tip 
(core) and propagate in radial direction. The aim of this 
paper is to investigate the effects of perturbing sources 
of nonlinear waves and their implications for the dynam- 
ics of spiral waves. We employ the complex Ginzburg- 
Landau equation (CGLE), which provides an universal 
description of spatially extended oscillatory systems near 
a supercritical Hopf bifurcation 01 . 

In this framework, we perturb a source of periodic 
waves in ID by moving its position back-and-forth in 
space. This motion of the source leads to a modulation 
in amplitude, wavelength and frequency of the emitted 
waves. We find that the modulation of the nonlinear 
waves is uniquely determined by the temporal period of 
the source motion in contrast to linear waves emitted by 
a moving source, where the source velocity is the rele- 
vant quantity. The richest scenario is found if the emit- 
ted waves are convectively Eckhaus unstable. If we con- 
sider a periodically moving source in the latter case, the 
modulations in amplitude and period of the waves may 
(a) get exponentially damped, (b) saturate or (c) grow 
monotonously far away from the source depending on the 
frequency of the applied forcing. 

Such a periodically moving source reproduces the ra- 
dial dynamics of rotating spiral waves subject to external 
forcing ^ or to the frequentlyobserved oscillatory insta- 
bility termed "meandering" In 2D, the radial mod- 
ulations caused by meandering lead to a second spiral 
superimposed on the simple rotating spiral; the resulting 
structure is called superspiral (see Fig. 1). Our study 
provides insight in the nonlinear behavior of such super- 
spirals extending earlier work based on a linear analy- 
sis . Unbounded growth of the modulation will lead to 
occurrence of space-time defects in ID and superspiral 
breakup in 2D, provided the system is sufficiently large. 




FIG. 1: Simple spiral waves (a-c) and superspirals (d-f) are 
observed in experiments (a,d) of the BZ system jg and in 
numerical simulations of the CGLE (b,e) show Re[A] 

and (c,f) |Aj. (b,c) /i = 1 and (e,f) inhomogeneous /i = 
l+0.7exp(— r^/10). Other parameters are ci = 3.5, ca = 0.34, 
system size is 512 x 512. 



Experimental evidence of such behavior was found in the 
Belousov-Zhabotinsky (BZ) reaction @. The first exper- 
imental observation of a superspiral was reported for a 
spiral subject to strong external forcing near its core |^. 

Complex Ginzburg-Landau equation. - Consider a spa- 
tially extended oscillatory medium described by the com- 
plex Ginzburg-Landau equation P, 0| 

dtA = fiA + (1 + ici)AA - (1 - ic3)\A\^A (1) 

with /i = 1. The complex field A{r,t) gives ampli- 
tude and phase of local oscillations depending on real 
coefficients ci , C3 determined by the underlying specific 
model at the onset of oscillations. Eq. Q exhibits 



plane wave solutions A{r, t) — -\/l — g2gi(9r-wt) 
u) = — C3+(7^(ci-|-C3) in an infinite or periodic ID medium 
with A = d^. For fixed control parameters, these waves 
become unstable if the wavenumber q is bigger than the 
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FIG. 2: Doppler effect of nonlinear waves due to periodic 
oscillations of the source at the left boundary. Eq. was 
integrated for increasing forcing periods (a) r = 8, (b) r = i3, 
(c) T = 15 with Eq. © and ci = 3.5, C3 = OA,Rs = 1. 
Dotted lines are guides to the eye that converge to |yl|min 
and \A\^ax in (b). 



Eckhaus wavenumber g|; — (1— ciC3)/(2(l-|-c§)-|-l— C1C3). 
The plane waves represent a one-parameter family that is 
parametrized by the wavenumber q. A source of periodic 
waves at r = is easily realized by applying the BCs 



A{r = 0) = and drA{r = L) = 0. 



(2) 



The specification of these BCs leads to the selection of a 
unique wavenumber qs- In this case, the exact solution 
of the CGLE is of the form 



A(r,t)=F(r)e^(/W-^*), 



(3) 



with the asymptotic behavior df{r)/dr — > qs,F{r) 
\/i — Qs for r — > 00 (far-field) and df{r)/dr ^ r, F{r) ^ 
r for r (near-field). An analytic expression for qs 
has been derived 0, 

= (ci -I- C3)g| -I- 3a{ci,C3)qs - C3 - 2ciq;(ci, 03)^, (4) 

where a is a function of the control parameters ci , C3 . 

In 2D, Eq. with A = 9^ + l/rdr + l/r'^dg possesses 
rotating spiral solutions with A = in the spiral core and 
in the far-field a selected wavenumber qs , which is similar 
to (4) as verified numerically ^, 10] . Eq. (4) allows to 
discriminate if the source or the spiral core emits stable 
(Eckhaus unstable) wavetrains with qs < qs (qs > <1e)- 
Further analysis showed that the Eckhaus instability is 
of convective nature and becomes absolute for sufficiently 
large values of ci and C3 |0. At given parameters, we 
may define a wavenumber qA > qs that characterizes 
absolutely unstable waves with q > qA- In finite systems, 
spirals were found to be stable as long as qs < qA 12,.13.] . 

Simulation results. - Here, we analyze the effects of an 
instability or of an external perturbation at the source 
or at the spiral core (meandering). We choose to per- 
turb the source in Eq. (j^J by varying its position rg. 
Consequently, 

^{f ^ fs) — with rs = Rs cos(27rt/T) 

and drA{r = L) = 0, (5) 



is used as EC of Eq. ^ in ID. The core motion of a me- 
andering spiral can be viewed as a source moving along 
a circle in 2D. A projection onto a radial direction for 
fixed angle yields a periodic back-and-forth motion. A 
sinusoidal motion like in Eq. (5) is expected near on- 
set, where the normal form of meandering [T^ provides 
a valid description. 

Hence, the BC in Eq. (5) also captures the radial dy- 
namics of a meandering spiral. In the numerical simu- 
lations the moving source initially modulates the local 
wavenumber which subsequently leads to a modulation 
of the wave amplitude. For large r (far-field), we observe 
modulations with a period T that is equal to the forc- 
ing period t and independent from Rs- First, we forced 
sources (respectively "spirals") in the parameter region 
where qs < qs, in this case the perturbation is always 
damped out while the waves move away from the source. 
A larger variety of responses occurs in the convectively 
Eckhaus unstable regime, where the unperturbed source 
selects a wavenumber qs in the interval [qE,qA]- Fig- 2 
shows three qualitatively different resulting profiles of the 
amplitude \ A\ as a function of the radial coordinate. Figs. 
2(a), (b) and (c) correspond to increasing values of the 
forcing period r. For small r-values, the result is similar 
to the case with qs < qs, the amplitude modulation is 
exponentially damped and is barely visible at sufficiently 
large r (Fig. 2(a)). For intermediate r, the amplitude 
modulation first grows and then reaches saturation (Fig. 
2(b)). The profile of \A\ is periodic and travels with a 
non-zero velocity; in the far-field, the radial dynamics 
resemble a so-called modulated amplitude wave (see be- 
low). Finally for large t, the amplitude modulation grows 
monotonically with r and space-time defects are formed 
at the breakup radius r = Rbu ~ 950 (Fig. 2(c)). For 
simulations with qs > qA, we always observe monotonous 
growth of the modulation leading to space-time defects. 

A drawback of the homogeneous CGLE in 2D is that it 
does not support meandering spirals. A non-saturating 
meandering instability has been observed for large values 
of ci [i5| . Alternatively, addition of a heterogeneity near 
the spiral core leads to meandering behavior similar to 
the one typically seen in reaction-diffusion systems [T^ . 
Figs. 1 (e),(f) show such a superspiral. The correspond- 
ing quantities for a regular spiral are shown in Fig. 1 
(b),(c). As in Fig. 2(b) and in the recent experiment by 
Zhou and Ouyang (see Fig. 1(d) and '^) the amplitude 
modulation shown in Fig. 1(f) saturates in the far-field. 
In the following, the study will be limited to the radial 
dynamics. One of the new results of our analysis is that 
the modulation of the amplitude may eventually satu- 
rate (see Fig. 2(b)). It also travels with constant shape 
and with a speed different from that of the underlying 
phase waves. To summarize, periodically moving sources 
emit modulated amplitude waves analogous to the way 
stationary sources send out plane waves. 

Modulated amplitude waves (MAWs) and superspirals. 
- MAWs are solutions of the CGLE and have the follow- 
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FIG. 3: (a) Bifurcation diagram for MAWs with period T and 
existence domains of spirals, superspirals and spiral breakup 
for various forcing periods r. Curves were computed as in 
IT^ and symbols denote minima and maxima of |j4| as 
measured in simulations, parameters are the same as in Fig. 
2. Solid (dashed) curves correspond to stable (unstable) solu- 
tions. Arrows (same for |A|maa:) indicate the evolution of ini- 
tial perturbations as they move away from the source. Their 
asymptotic values depend on T = r but not on Rs- (b) 
Breakup radius Kbu at which defects firstly occur in sim- 
ulations. Smaller remaining spirals result from larger Ks 
{Rs = 5 circles, -Rs = 1 squares). 



ing general form 

A{r, t) 



{z)e 



ujt) 



(6) 



with unknown periodic functions a and <f> of the comov- 
ing coordinate z = r — vt 17J. The analysis of MAWs 
has revealed that they bifurcate from plane waves with 
q > qe- MAWs form a two-parameter family of waves 
described by the wavelength P of the modulation of the 
amplitude a = |j4| and by the spatially averaged phase 
gradient u — q + {dr(p). The velocity v and the period 
T = P/v can be derived from these quantities. The aver- 
age phase gradient is conserved as long as no space-time 
defect is formed. For periodically perturbed sources, the 
phase gradient value is fixed to v = qs- Therefore, only 
one free parameter is left. For the present purpose it is 
most convenient to use the period T. The resulting bi- 
furcation diagram is shown in Fig. 3 for the case qs > Qe 
(Eckhaus unstable range). In the Eckhaus stable regime, 
no stable MAWs exist. 

Fig. 3 reveals that stable MAWs exist for periods T 
in the interval [ThbtTsn]- They are "born" in a Hopf 
bifurcation (HB) and "die" in a saddle-node bifurcation 
(SN). For T < Thb, the plane wave is stable and for 
T > TsN no MAWs exists and the dynamics leads to the 
formation of defects 0,^3- This clarifies our findings in 
Fig. 2: the profiles in Fig. 2(a), 2(b) and 2(c) correspond 
to forcing within the regions r < Thb, Thb < t < Tsn 
and Tsn < t, respectively. Fig. 3 (b) shows the breakup 
radius Rbu found in simulations for periods t > Tsn- 
It is crucial to realize that the properties of the sat- 




FIG. 4: Phase diagram indicating the region (shaded) where 
superspirals may occur in the CGLE for "rocking" sources. 
Simple spirals with selected qs (01,03) are convectively un- 
stable between the thick solid curves. Superspiral breakup 
is only possible between the dashed and dash-dotted curves. 
The Eckhaus instability for arbitrary q is supercritical above 
the thin solid curve and the thin dotted curve indicates the 
Benjamin- Feir-Newell line (BFN). 



urated modulations, observed in numerical simulations 
with a periodically moving source, indeed correspond to 
the MAWs with v = qs and T = t computed via a bi- 
furcation analysis (as shown in Fig. 3). Thus, a unique 
MAW characterized by two parameters is selected in ID 
by BC in Eq. (5) or in 2D by the two intrinsic frequen- 
cies of a meandering spiral. We conjecture, that the bi- 
furcation diagram of superspirals in 2D can be predicted 
from the corresponding bifurcation diagram of the far- 
field MAWs provided the meandering period (that here 
corresponds to r) is known. Hence, the superspirals with 
saturated modulation should cease to exist in the saddle- 
node bifurcation of the associated MAW. If r > Tsn, 
the modulations grow monotonously and lead to appear- 
ance of space-time defects in ID, respectively topological 
defects in 2D (spiral breakup). In explicit 2D simula- 
tions, such behavior has already been observed in a ho- 
mogeneous reaction-diffusion model for calcium waves, 
where simultaneous appearance of the Eckhaus instabil- 
ity and meandering leads to a breakup far away from 
the spiral core [23|. Future studies may address the pos- 
sibility of MAWs in realistic reaction-diffusion systems 
which exhibit Eckhaus instabilities ^0, 21]. A previous 
linear analysis of superspirals predicted that the mean- 
dering instability of a spiral with an Eckhaus unstable 
wavetrain in the far-field may produce superspirals with 
exponentially growing modulation, while in the standard 
case of a meandering spiral emitting a stable wave train 
exponential damping of the modulation should be ob- 
served Our nonlinear analysis introduces superspi- 
rals with a saturated modulation in the far-field as a 
third possibility, thus supporting the experimental ob- 
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servation of such structures in the BZ reaction (Fig. 1 
and 0). We suggest an interpretation of the reported 
scenario damped superspiral - saturated superspiral - far- 
field breakup: the three phenomenologies simply corre- 
spond to the three mentioned regions of the bifurcation 
diagram for the MAWs. Notice that bifurcation diagrams 
similar to Fig . 3 (a) are found when ci or C3 are varied 
at fixed T [ll 0,113 • 

Phase diagram. - Finally we determine the regions 
of the ci — C3 phase diagram where the discussed 
phenomena may be found. The results are summarized 
in Fig. 4 where v is fixed equal to qs {01,03) as given 
by Eq. (4). Saturated superspirals can appear only 
above the Eckhaus line for qs > qs- For small values 
of ci, the MAWs bifurcate subcritically and are always 
unstable Fig. 4 also shows the line of absolute 

instability [qs — qA) above which simple spirals break 
up giving rise to spatio-temporal chaos. Stable MAWs 
with v — qs do exist in the shaded re gion (for an ex- 
tensive discussion on MAW stability see [13 )■ Below the 
thick dashed line in Fig. 4 no saddle-node bifurcation 
does occur and breakup is therefore prevented. Between 
the thick dashed and the thick dash-dotted lines in Fig. 
4, bifurcation diagrams similar to Fig. 3(a) are found. 
In this region, all the three behaviors reported in Fig. 
2 are possible, depending on the period of the forcing 
or meandering. Altogether, the results here presented 
link far-field breakup of meandering spirals to a saddle- 



node bifurcation. Thus, this route to spatio-temporal 
chaos can be distinguished from the previously reported 
scenario of far-field breakup, where spiral breakup has 
been explained by an absolute Eckhaus instability of the 
asymptotic far-field wavetrain [13. ITsL f20l |23| . 

Conclusions. - We have studied the Doppler effect 
associated with the back-and-forth motion of a source 
emitting periodic nonlinear waves. Usually, the result- 
ing modulation dies out by exponential damping as the 
waves move away from the source. If the emitted wave- 
train is convectively Eckhaus unstable, the wave modu- 
lation can also saturate or grow exponentially far away 
from the source depending on the period of the back- 
and-forth motion. These scenarios are fully determined 
by the bifurcation diagram of corresponding MAWs. Our 
results offer a consistent explanation of recent experimen- 
tal results obtained in a chemical reaction S]. Moreover, 
they may find application in future studies of other ex- 
perimental systems exhibiting convective instabilities like 
hydrothermal waves . Periodic forcing of sources near 
the transition to chaotic or turbulent dynamics may be 
used to probe the existence and properties of modulated 
structures. 
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